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A search for the rare decays Bg — ?► /^+/^ and B^ — )■ /^+^ is performed in pp collisions 
at ^/s = 7 TeV, with a data sample corresponding to an integrated luminosity of 5 fb^ 
collected by the CMS experiment at the LHC. In both decays, the number of events 
observed after all selection requirements is consistent with the expectation from back- 
ground plus standard model signal predictions. The resulting upper limits on the 
branching fractions are e(B° -^ }i^}i') < 7.7 x 10"'^ and i3(B0 ^ fi+ji-) < 1.8 x 10"'' 
at 95% confidence level. 
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1 Introduction 

The decays Bs(B°) — > fi^f/^ are highly suppressed in the standard model (SM) of particle 
physics, which predicts the branching fractions to be S(Bg — > }i^}i^) = (3.2 ± 0.2) x 10^^ 
and /5(B° — > fi^f/^) = (1.0 ± 0.1) x 10^^° |lj. This suppression is due to the flavor-changing 
neutral current transitions b — ?■ s(d), which are forbidden at tree level and can only proceed via 
high-order diagrams that are described by electroweak penguin and box diagrams at the one- 
loop level. Additionally, the decays are helicity suppressed by a factor of m^^/m^, where m^, 
and niQ are the masses of the muon and B meson, respectively (the symbol B is used to denote 
B*^ or Bg mesons). Furthermore, these decays also require an internal quark annihilation within 
the B meson that reduces the decay rate by an additional factor of /^ /m|, where /b is the decay 
constant of the B meson. The leading theoretical uncertainty is due to incomplete knowledge 
of /b, which is constrained by measurements of the mixing mass difference Amg (Am^) for Bg 
(B*^) mesons. 

Several extensions of the SM predict enhancements to the branching fractions for these rare 
decays. In supersymmetric models w^ith non-universal Higgs masses [2J and in specific models 
containing leptoquarks ||3|, for example, the Bg — )■ }i'^}i^ and B*^ — > }i^}i^ branching fractions 
can be enhanced. In the minimal supersymmetric extension of the SM, the rates are strongly 
enhanced at large values of tan /3, w^hich is the ratio of the two vacuum expectation values of 
the two Higgs boson doublets ||ll|5j . However, in most models of new physics, the decay rates 
can also be suppressed for specific choices of model parameters ||6l. 

At the Tevatron, the DO experiment has published an upper limit of B(Bg — > ji^jx^) < 5.1 x 
10"^ [7J at 95% confidence level (CL). The CDF experiment has set a limit of 6(B° -^ fi^H') < 
4.0 X 10^^ and 6(B° —> }i^}i^) < 6.0 x 10^^, and also reported an excess of Bg —> }i^}i^ events, 
corresponding to B{Bl -^ ^^^') = (l-S+Jg) x 10"^ |8|. At the Large Hadron Collider (LHC), 
two experiments have published results: S(Bg -^ H^H^) < 1-9 x 10^^ and i5(B° -^ H^h) < 
4.6 X lO"'^ by the Compact Muon Solenoid (CMS) Collaboration [9], and B(B° -^ ^^}l') < 
1.4 X 10-s and B{B° -^ ^^^') < 3.2 x lO"'' by the LHCb Collaboration IflOH. 

This paper reports on a new simultaneous search for Bg — > }i^}i' and B*^ — > }i'^}i^ decays us- 
ing data collected in 2011 by the CMS experiment in pp collisions at v^ = 7TeV at the LHC. 
The dataset corresponds to an integrated luminosity of 5 fb^ . An event-counting experiment 
is performed in dimuon mass regions around the Bg and B*^ masses. To avoid potential bias, 
a "blind" analysis approach is applied where the signal region is not observed until all selec- 
tion criteria are established. Monte Carlo (MC) simulations are used to estimate backgrounds 
due to B decays. Combinatorial backgrounds are evaluated from the data in dimuon invariant 
mass {m^m) sidebands. In the CMS detector, the mass resolution, which influences the sepa- 
ration between Bg — > }i^}i^ and B° — >• }i^}i~ decays, depends on the pseudorapidity rj of the 
reconstructed particles. The pseudorapidity is defined as t] = — ln[tan(0/2)], where is the 
polar angle with respect to the counterclockwse proton beam direction. The background level 
also depends significantly on the tj of the B candidate. Therefore, the analysis is performed 
separately in two channels, "barrel" and "endcap", and then combined for the final result. The 
barrel channel contains the candidates where both muons have \r]\ < 1.4 and the endcap chan- 
nel contains those where at least one muon has \f]\ > 1.4. 

A "normalization" sample of events with B+ — )■ J/j/^K^ decays (where J/i/? — ?► fi^ji^) is used to 
remove uncertainties related to the bb production cross section and the integrated luminosity. 
The signal and normalization efficiencies are determined through MC simulation studies. To 
validate the simulation distributions, such as the Bg transverse momentum (pj) spectrum, and 
to evaluate potential effects resulting from differences in the fragmentation of B+ and Bg, a 
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"control" sample of reconstructed Bg —> ]/ip(p decays (with ]/ip —> }l^}i^ and (p —> K+K^) is 
used. 

The dataset includes periods of high instantaneous luminosity conditions, with an average of 8 
interactions per bunch crossing (later referred to as "pileup"). The analysis algorithms and the 
selection criteria have been optimized to mitigate the effects of pileup by reducing the influence 
of tracks coming from additional interactions in the event, as explained in SectionlSl In parallel 
with the LHC luminosity increase, the CMS event triggering requirements also changed during 
the data-taking period. The analysis and simulations take these changes into account so that 
all MC samples incorporate the appropriate mixture of the trigger conditions, and the selection 
requirements applied in the data reconstruction are more restrictive than the most stringent 
trigger criteria. 

The limits on the branching fractions depend on both systematic and statistical uncertainties. 
Several sources of systematic uncertainties can influence the estimated efficiency: detector ac- 
ceptance, and analysis, muon identification and triggering efficiencies. The evaluation of the 
individual values are presented in the sections below when discussing the relevant efficiencies 
and then are combined in Section |6l 

The data analyzed here include the event sample corresponding to an integrated luminosity of 
1.14 fb^ , w^hich w^as used to obtain the earlier CMS result ||9l. The present analysis differs in 
several ways: the total dataset is almost five times larger; new^ selection variables are added 
to the analysis; the selection criteria are optimized for higher pileup and varying trigger re- 
quirements; and the description of rare backgrounds is improved. All these changes result in a 
better signal sensitivity. 

2 Monte Carlo simulation 

Simulated events are used to determine the efficiencies for the signal and normalization sam- 
ples. We split the efficiency into four parts: detector acceptance, analysis efficiency, and muon 
identification and trigger efficiencies. The detector acceptance combines the geometrical detec- 
tor acceptance and the tracking efficiency, and is defined as tracks within \i]\ < 2.4 and satisfy- 
ing pr > 1 GeV (pt > 0.5 GeV) for muons (kaons). The acceptance is about 25% (23%) for signal 
events in the barrel (endcap) channels. In the pj range relevant for this analysis the tracking 
efficiency for isolated muons and kaons is above 99.5% [llj. The analysis efficiency refers to 
the selection requirements described in Section |5} and is for signal events about 2.0% (1.2%) 
in the barrel (endcap) channels. The muon identification and trigger efficiencies are presented 
in Sections [3] and [4l respectively. The analysis, muon identification, and trigger efficiencies are 
all obtained from simulation and checked in data. Good agreement is found, and the residual 
differences are used to estimate systematic uncertainties on the efficiency estimates. 

The simulated samples are also used to estimate the background from rare B decays where one 
or two hadrons are misidentified as muons. These decays include a variety of channels of the 
type B — > h^/^+v and B — > h+h^, w^here h is a tt, K or p and B stands for B'^,Bg mesons or A\, 
baryons. The most important backgrounds are from B^ — > K^K+, B*^ — > K+tt^ and from the 
semileptonic decays B*^ — > iz^ji^v, Bg -^ K^/^+v, and A^ -^ pji^v. 

The samples of simulated events are generated with Pythia 6.424 (TUNE Z2) fT2], the un- 
stable particles are decayed via EvtGen IIT3l , and the detector response is simulated with 
Geant4 [14]. The signal and background events are selected from generic quantum chromo- 
dynamic (QCD) 2 — > 2 sub-processes and provide a mixture of gluon-fusion, flavor-excitation. 



and gluon-splitting production. The evolution of the triggers used to collect the data is incor- 
porated in the reconstruction of the simulated events. The number of simulated events in all 
the channels approximately match the expected number given the integrated luminosity. 

3 The CMS detector 

The CMS detector is a general-purpose detector designed and built to study physics at the 
TeV scale. A detailed description can be found in Ref. |15|. For this analysis, the main sub- 
detectors used are a silicon tracker, composed of pixel and strip detectors within a 3.8 T axial 
magnetic field, and a muon detector, which is divided into a barrel section and two endcaps, 
consisting of gas-ionization detectors embedded in the steel return yoke of the solenoid. The 
silicon tracker detects charged particles within the pseudorapidity range \r]\ < 2.5. The pixel 
detector is composed of three layers in the barrel region and two disks located on each side 
in the forward regions of the detector. In total, the pixel detector contains about 66 million 
100 }im X 150 jim pixels. Further from the interaction region is a microstrip detector, which is 
composed of ten barrel layers, and three inner and nine outer disks on either end of the detector, 
with a strip pitches between 80 and 180 pim. In total, the microstrip detector contains around 
10 million strips and, together with the pixel detector, provides an impact parameter resolution 
of ~ 15 f/m. Due to the high granularity of the silicon tracker and to the strong magnetic field, 
a Pt resolution of about 1.5% [16] is obtained for the charged particles in the pr range relevant 
for this analysis. The systematic uncertainty on the hadronic track reconstruction efficiency is 
estimated to be 4% ||T6|. Muons are detected in the pseudorapidity range \ri\ < 2.4 by detectors 
made of three technologies: drift tubes, cathode strip chambers, and resistive plate chambers. 
The analysis is nearly independent of pileup because of the high granularity of the CMS silicon 
tracker and the excellent three-dimensional (3D) hit resolution of the pixel detector. 

The dimuon candidate events are selected with a two-level trigger system, the first level only 
uses the muon detector information, while the high-level trigger (HLT) uses additional infor- 
mation from the pixel and strip detectors. The first-level trigger requires two muon candidates 
without any explicit pj requirement, but there is an implicit selection since muons must reach 
the muon detectors (about 3.5 GeV in the barrel and 2GeV in the endcap). The HLT imposes a 
Pt requirement and uses additional information from the silicon tracker. As the LHC instan- 
taneous luminosity increased, the trigger requirements w^ere gradually tightened. This change 
in trigger requirements is also included in the trigger simulations. The most stringent HLT 
selection requires two muons each with pj > 4 GeV, the dimuon pj > 3.9 GeV (5.9 GeV in the 
endcap), dimuon invariant mass within 4.8 < m^j^^ < 6.0 GeV, and a 3D distance of closest 
approach to each other of d'^^ < 0.5 cm. For the entire dataset, the offline analysis selection is 
more restrictive than the most stringent trigger selections. 

For the normalization (B+ — > J/t/;K+) and control (Bg — > ]/ip(p) samples, the data are collected 
by requiring the following: two muons each with pj > A GeV, dimuon pj > 6.9 GeV, \}]\ < 2.2, 
invariant mass within 2.9 < m^ifi < 3.3 GeV, d^^ < 0.5 cm, and the probability of the x^ per 
degree of freedom (;)|;^/dof) of the dimuon vertex fit greater than 15%. To reduce the rate of 
prompt ]/ip candidates, two additional requirements are imposed in the transverse plane: (i) the 
pointing angle Uxy betw^een the dimuon momentum and the vector from the beamspot (defined 
as the average interaction point) to the dimuon vertex must fulfill cos ocxif > 0.9; and (ii) the 
flight distance significance ixy/ci^xy) niust be larger than 3, where £xy is the two-dimensional 
distance between the primary and dimuon vertices and cr{£xy) is its uncertainty. 

The trigger efficiencies for the various samples are determined from the MC simulation. They 
are calculated after all muon identification selection criteria, as discussed in Section |4} have 
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been applied. For the signal events the average trigger efficiency is 84% (74%) in the barrel 
(endcap) channel. The trigger efficiency for the normalization and control samples varies from 
77% in the barrel channel to 60% in the endcap channel. This analysis depends on the ratio 
of the signal efficiency to the normalization sample efficiency. The systematic uncertainty on 
the trigger efficiency ratio is estimated as the sum in quadrature of two components. The 
first component is defined as the variation of the efficiency ratio when varying the muon p^ 
threshold from 4 to 8 GeV in the MC simulation. The second one is the difference between the 
ratios determined in data and MC simulations using the tag-and-probe approach (described in 
Section [4]|. The systematic uncertainty on the ratio is estimated to be 3% in the barrel channel 
and 6% in the endcap channel. 

4 Muon identification 

Muon candidates are reconstructed by combining tracks found in the silicon tracker and the 
muon detector lllTl IT8| . In order to ensure high-purity muons, the following additional re- 
quirements are applied: (i) muon candidates must have at least two track segments in the 
muon stations; (ii) they must have more than 10 hits in the silicon tracker, of which at least 
one must be in the pixel detector; (iii) the combined track must have ,^^/dof < 10; and (iv) the 
impact parameter in the transverse plane d^y, calculated with respect to the beamspot, must be 
smaller than 0.2 cm. The systematic uncertainty on the muon track reconstruction efficiency is 
2% [TTl and is included in the uncertainty of the total efficiency. 

The ratio of the muon identification efficiencies between the signal and normalization samples 
is used in this analysis. This ratio is determined in tw^o ways. First, the MC event samples 
contain a full simulation of the muon detector, which allows an efficiency determination by 
counting the events that pass or fail the muon identification algorithm. Second, the muon 
identification efficiency is determined with a tag-and-probe method [17J, w^hich is applied to 
both data and MC event samples. To study the single-muon identification efficiency, the decays 
J/j/' — >• ]i^}i^ are used. In the tag-and-probe method, a "tag" muon, satisfying strict muon crite- 
ria, is paired with a"probe" track, where together they combine to give the J/t/^ invariant mass, 
thus indicating the probe is in fact a muon. The single-muon efficiency is determined by the 
number of probe tracks passing or failing the muon identification algorithm. Dedicated trigger 
paths constructed using the tag muon and either a silicon track or a signal in the muon cham- 
bers are employed for this study, which ensures large event samples while avoiding potential 
bias of the efficiency measurement from using events triggered by the probe. 

The muon identification efficiency is calculated after all selection criteria, including the detector 
acceptance, have been applied. For the signal events, the average efficiency is 71% (85%) in 
the barrel (endcap) channel based on the MC simulation. For the normalization and control 
samples, the muon identification efficiency is about 77% (78%) in the barrel (endcap). Pair- 
correlation effects influence these numbers |17|. The dimuon efficiency can be altered with 
respect to the product of single-muon efficiencies depending on the mutual proximity of the 
two muons in the muon system. This effect is included in the efficiency calculations in the 
detailed MC simulation of the muon detectors. The systematic uncertainty on the identification 
efficiency ratio is estimated in the same way as for the muon trigger efficiency ratio (SectionlSll, 
and is 4% in the barrel and 8% in the endcap. 



5 Analysis 

The reconstruction of B — > }i^}i^ candidates requires two oppositely-charged muons that orig- 
inate from a common vertex and have an invariant mass in the range 4.9 < m^^ < 5.9 GeV. A 
fit of the B-candidate vertex is performed and its x^/'ioi is evaluated. The two daughter muon 
tracks are combined to form the B-candidate track. 

The primary vertex associated with a B candidate is chosen from all reconstructed primary 
vertices as the one which has minimal separation along the z axis from the z intercept of the 
extrapolated B candidate track. Reconstruction effects due to pileup are largely eliminated 
by the primary vertex matching procedure. The position of this primary vertex is then refit 
without the tracks of the B candidate with an adaptive vertex fit fW\, where tracks are assigned 
a weight < w < 1 based on their proximity to the primary vertex. After the refit, B candidates 
•with badly reconstructed primary vertices are eliminated by requiring the average track weight 
{w) > 0.6. The 3D impact parameter of the B candidate S^d, its uncertainty ^{S^u), and its 
significance ^3d/c^(^3d) are measured with respect to the primary vertex. 

The isolation of the B candidate is an important criterion in separating the signal from back- 
ground. Three variables are used for this purpose: 

• The isolation variable I = Pt{B)/ {pt{B) + I^trkPr) is calculated from the trans- 
verse momentum of the B candidate pr(B) and the transverse momenta of all other 
charged tracks satisfying AR = y^{Arj)'^ + (A^)^ < 0.7, where At] and Ai^ are the 
differences in pseudorapidity and azimuthal angle between a charged track and the 
B-candidate momentum. The sum includes all tracks with pj > 0.9 GeV that are 
(i) consistent with originating from the same primary vertex as the B candidate or 
(ii) have a distance of closest approach dca < 0.05 cm with respect to the B vertex 
and are not associated w^ith any other primary vertex . 

• The number of tracks N^^^^^ with pj > 0.5 GeV and dca < 0.03 cm with respect to the 
B-candidate's vertex. 

• The minimum distance of closest approach betw^een tracks and the B-candidate's 
vertex, d^^, for all tracks in the event that are either associated with the same primary 
vertex as the B-candidate or not associated with any other primary vertex. 

The first variable describes the isolation primarily with respect to tracks coming from the pri- 
mary vertex itself. The latter two variables quantify the isolation of the B vertex. They help to 
reject partly reconstructed B decays where there are other tracks in addition to the two muons 
associated with the B-candidate vertex. 

The distributions of the variables described above are shown in Fig. IT] for signal events from 
the MC simulation and for data background events. These include the momenta of the higher- 
momentum (leading) and lower-momentum (sub-leading) muons pr,;*! and pT,p2, Pt(B), the 
3D pointing angle a^^, the 3D flight length significance ^3d/c^(^3d)/ the ;\;^/dof, and the iso- 
lation variables {I, N^^^^^, and d^^). The data background events are defined as B candidates 
with a dimuon mass in the sidebands covering the range 4.9 < m^,^, < 5.9 GeV, excluding the 
(blinded) signal window^s from 5.20 < m^,^ < 5.45 GeV. Events shown in Fig. Ill must pass a 
tight selection that is close to the final one: muon pj > 4 GeV, Pt(B) > 7.5 GeV, a < 0.05, 
;^Vdof < 2, ^3d/o-(^3d) > 15, Ssd < 0.008cm, ^3d/o-('53d) < 2, J > 0.8, d^^ > 0.015cm, and 
^trk^*^ < 2 tracks. For each distribution, the selection requirements for all variables, apart from 
the one plotted, are applied. This figure illustrates the differences in the distributions of signal 
and background events, and shows which variables are effective in reducing the background 
events, e.g., ^3d/c^(^3d)- The analysis efficiency for each selection requirement is determined 
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from the simulated events. 

The reconstruction of B+ — > J/)/?K+ — > /^+f/^K+ and Bg — >• J/j/^^ — ?► /i+f/^K+K^ events is very 
similar to the reconstruction of B — )■ /^^fi^ events. Candidates with two oppositely-charged 
muons sharing a common vertex and with invariant mass in the range 3.0-3.2 GeV are recon- 
structed. The selected candidates must have a dimuon pj > 7 GeV. Then they are combined 
w^ith one or two tracks each assumed to be a kaon, with pj > 0.5 GeV and \f]\ < 2.4. The 3D 
distance of closest approach between all pairs among the three (four) tracks is required to be 
less than 0.1 cm. For Bg — > ]/tp(p candidates, the two assumed kaon tracks must form an invari- 
ant mass in the range 0.995-1.045 GeV and have AR(K+,K") < 0.25. The three (four) tracks 
from the decay are used in the vertex fit. All requirements listed above for B — > }i^}i^ events 
are also applied here, including the vertex-fit selection ;\;^/dof < 2, which eliminates poorly 
reconstructed candidates. Only B candidates with an invariant mass in the range 4.8-6.0 GeV 
are considered. 

Figures[2]and[3]show the MC simulation and sideband-subtracted data distributions for a num- 
ber of variables for the B+ — > J/j/?K+ and Bg — > ]/^(p candidates, respectively. For each distribu- 
tion, the selection requirements for all variables, apart from the one plotted, are applied. The 
relative efficiency for each selection requirement is determined separately in data and MC sim- 
ulation and compared. The largest relative differences are 2.5% for the isolation selection in the 
normalization sample and 1.6% for the ;\;^/dof selection in the control sample. We combine in 
quadrature the differences for all distributions to estimate the systematic uncertainty related to 
the selection efficiency and obtain 4% (3%) for the normalization (control) sample. The control 
sample uncertainty is used for the signal sample. 

The dataset used in this analysis is affected by pileup, which includes an average of 8 recon- 
structed primary vertices per event. The distribution of the primary vertex z position has a 
Gaussian shape with an RMS of approximately 5.6 cm. To study a possible dependence on 
the amount of pileup, the efficiencies of all selection criteria are calculated as a function of the 
number of reconstructed primary vertices. In Fig. [4] this dependence is shown for the normal- 
ization sample. A horizontal line is superimposed to guide the eye indicating that no significant 
dependence is observed. The same conclusion is also obtained in the MC simulation by com- 
paring the selection efficiency for events with less than six primary vertices to those with more 
than ten primary vertices. Similar studies of the control sample, albeit with less precision, lead 
to the same conclusion: the analysis is not affected by pileup. 

Variables sensitive to the underlying production processes (gluon fusion, flavor excitation, or 
gluon splitting) are also studied to validate the production process mixture in the MC simu- 
lation. The clearest distinction among the three processes is obtained by studying (i) the AR 
distribution between the B candidate and another muon and (ii) the pj spectrum of this muon. 
The MC simulation (Pythia) describes these distributions adequately. 

6 Results 

The present analysis differs significantly from the previous one HI: 

• The muon identification algorithm has changed. A tighter selection is used, which 
significantly decreases the rate at which kaons and pions are misidentified as muons. 

• The definition of isolation is different and tw^o additional isolation variables are 
used, w^hich reduces the influence of event pileup and also lowers the combinato- 
rial background. 
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Figure 1: Comparison of simulated Bg — )■ }i^}i^ decays and background dimuon distributions 
as measured in the mass sidebands. Top row: the transverse momentum for the leading muon, 
sub-leading muon, and B-candidate; middle row: the 3D pointing angle, flight length signifi- 
cance, and B-candidate's vertex ;ii;^/dof; bottom row: the isolation variables I, N^l^^, and d^^. 
The MC histograms are normalized to the number of events in the data. 
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Figure 2: Comparison of measured and simulated B^^ — ?► J/i/;K+ distributions. Top row: the 
transverse momentum for the leading muon, sub-leading muon, and B-candidate; middle row: 
the 3D pointing angle, flight length significance, and B-candidate's vertex ;^^/dof; bottom row: 
the isolation variables I, N^^^^*^, and d'^^. The MC histograms are normalized to the number of 
events in the data. 
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Figure 3: Comparison of measured and simulated Bg — ?► J/i/'^ distributions. Top row: the 
transverse momentum for the leading muon, sub-leading muon, and B-candidate; middle row: 
the 3D pointing angle, flight length significance, and B-candidate's vertex ;^^/dof; bottom row: 
the isolation variables I, N^^^^*^, and d'^^. The MC histograms are normalized to the number of 
events in the data. 
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• The requirement, for the normalization and control samples, that the two muons 
bend away from each other is removed, making the selection of these samples more 
similar to that for the signal. 

• The rare backgrounds, discussed below, are taken into account when calculating the 
combinatorial background, thus improving the background estimate in the signal 
window. 

The variables discussed in SectionlSlare optimized to obtain the best expected upper limit using 
MC signal events and data sideband events for the background. The optimization procedure 
is based on a random-grid search of about 1.4 x 10^ analysis selections. During this search, 
eleven variables are randomly sampled w^ithin predefined ranges. The resulting optimized 
requirements, which are used to obtain the final results, are summarized in Table [l] These 
requirements were established before observing the number of data events in the signal region. 
Hence, the analysis was blind to the signal events in the 5.20 < m^^, < 5.45 GeV mass range. In 
the endcap regions the selection is tighter than in the barrel because of the substantially larger 
background. The signal efficiencies £tot for these selections are shown in Table [2] They include 
all selection requirements: the detector acceptance, and the analysis, muon identification, and 
trigger efficiencies. The quoted errors include all the systematic uncertainties. In general, the 
present analysis uses more strict selection requirements than in the earlier analysis |9|, resulting 
in a higher sensitivity and a better signal-to-background ratio, but also a lower signal efficiency. 
As an additional test, the optimization was repeated to maximize the ratio S / \/S + B, where 
S is the number signal events and B is the number of background events. This resulted in 
a similar set of parameters to the ones listed in Table [ll but without an improvement in the 
expected upper limit. 

To evaluate a possible bias due to the optimization of the selection criteria in the data side- 
bands and to validate our background expectation, the following crosscheck is performed. All 
candidates with I < 0.7, including those within the blinded region, are selected ("inverted iso- 
lation" selection), which generates a background-enriched sample w^ith a very small expected 
signal contribution. From this sample, the candidate yields in the sidebands and in the blinded 
region are determined. The sideband yields are used to predict, through interpolation, the 
number of background candidates in the blinded region. Then the number of predicted back- 
ground events can be compared to the number of observed candidates in the blinded region. 
This comparison is performed separately for the barrel and endcap channels and for the Bg and 
B'^ signal windows. Within statistical uncertainties, good agreement is found for all four cases, 
which means that no significant biases are present in the background interpolation. 

The simulated dimuon mass resolution for signal events depends on the pseudorapidity of the 
B candidate and ranges from 37MeV for ?/ ~ to 77MeV for |//| > 1.8. The dimuon mass scale 
and resolution in the MC simulation are compared with the measured detector performance 
by studying ]/xp — > }i^}i^ and Y(1S) — > }i^}i^ decays. The residual differences between sim- 
ulation and data are small and the uncertainty on the efficiency coming from these effects is 
estimated to be 3%. 

Branching fractions are measured separately in the barrel and endcap channels using the fol- 
lowing equation 

S(B°^^V) = ^4^^S(B+)^ (1) 

Nobs /s etot 

where etot is the total signal efficiency, ^^^^ is the number of reconstructed B+ — > J/t/jK+ decays, 
e^ot is the total efficiency of B+ reconstruction, i3(B+) is the branching fraction for B+ — > 
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6 Results 



Table 1: Selection criteria for B 







^+^- 



and B^ — > fi^fi search in the barrel and endcap. 



Variable 


Barrel 


Endcap 


units 


PT,a > 


4.5 


4.5 


GeV 


VT^,l > 


4.0 


4.2 


GeV 


Ptb> 


6.5 


8.5 


GeV 


hu < 


0.008 


0.008 


cm 


^3d/^(^3d) < 


2.000 


2.000 




DC < 


0.050 


0.030 


rad 


X^/doi < 


2.2 


1.8 




i3d/(r{i3d) > 


13.0 


15.0 




I> 


0.80 


0.80 




4.> 


0.015 


0.015 


cm 


^close < 


2 


2 


tracks 



J/t/'K^ 



/^+/^ K 



/ /u//s is the ratio of the B+ and Bg production cross sections, and Ns is 



the background-subtracted number of observed Bg —> fi^fi^ candidates in the signal window 
5.30 < nifipi < 5.45 GeV. The width of the signal windows is adjusted to maximize the ef- 
ficiency for the Bg — > f/+fi^ decay, and it is approximately equal to twice the expected mass 
resolution in the endcap region. We use the value /s//u = 0.267 ± 0.021, measured by LHCb 
for 2 < ?/ < 5 [19J and 5(B+) = i3(B+ -^ J/1/7K+ -^ /^+fi"K+) = (6.0 ± 0.2) x 10"^ [20J. 
An analogous equation is used to measure the B^ — ?► }i^}i^ branching fraction, w^ith the signal 
window 5.2 < m^^ < 5.3 GeV and the ratio /d//u taken to be 1. 

The number of reconstructed B+ mesons N^^,^ is (82.7 ± 4.2) x lO'^ in the barrel and (23.8 ± 
1.2) X 10'^ in the endcap. The invariant mass distributions are fit w^ith a double-Gaussian func- 
tion for the signal and an exponential plus an error function for the background, as shown 
in Fig. pi Partially reconstructed B^ decays (e.g., B — )■ ]/tpK* with one of the K* decay prod- 
ucts not reconstructed) lead to a step function-like behavior at a mass of m ?a 5.15 GeV. This 
background shape was studied in detail in MC simulation and is parametrized with an error 
function of different width in the barrel and endcap. The systematic uncertainty on the fit 
yield, 5% in the barrel and in the endcap, is estimated by considering alternative fitting func- 
tions and by performing a fit with the dimuon invariant mass constrained to the J/t/; mass. The 
total efficiency ef^j, including the detector acceptance, is determined from MC simulation to 
be (11.0 ± 0.9) X 10"^ for the barrel and (3.2 ± 0.4) x 10"* for the endcap, where the errors 
include statistical and systematical uncertainties. The detector acceptance part (which includes 
the track finding efficiency) of the total efficiency has a systematic uncertainty of 3.5% (5.0%) 
in the barrel (endcap). It is estimated by comparing the values obtained separately w^ith three 
different bb production mechanisms: gluon splitting, flavor excitation, and flavor creation. 

The branching fraction for the control decay Bg — > ]/tp(p which was analyzed in parallel with 
the normalization and signal decays, has also been evaluated using Equation IT] The resulting 
branching ratio is in agreement with the world average [20] . Moreover, the results for the 
barrel and endcap channels agree within the statistical uncertainties, showing the validity of 
extending the /s //u measurement from |IT9| to the barrel region. 



Events in the signal w^indow have several sources: (i) genuine signal decays, (ii) decays of the 
type B — > hh', where h, h' are charged hadrons misidentified as muons (referred to as "peaking" 
background), (iii) rare semileptonic decays B — > hjiv, where h is misidentified as a muon, 
and (iv) combinatorial background. Note that events from the third category predominantly 
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populate the lower sideband. 

The expected numbers of signal events N^g^ai ^°^ *^^ barrel and endcap channels are shown in 
Table |2] They are calculated assuming the SM branching fractions lU and are normalized to 
the measured B+ yield. 

The expected numbers of rare semileptonic decays and peaking background events, N^^^, are 
also shown in Table [2] They are evaluated from a MC simulation, which is normalized to the 
measured B+ yields, and from muon misidentification rates measured in D*+ — > D^n^, D^ — > 
K^ 7T+, and A — > p/r" samples IITtII . The average misidentification probabilities in the kinematic 
range of this analysis are (0.10 ± 0.02)% for pions and kaons, and (0.05 it 0.01)% for protons, 
where the uncertainties are statistical. The systematic uncertainty on the background includes 
the uncertainties on the production ratio (for Bg and A^ decays), the branching fraction, and 
the misidentification probability. 

Also shown in Table 2 are the expected numbers of combinatorial background events N^q^^- 
They are evaluated by interpolating into the signal window^ the number of events observed 
in the sideband regions, after subtracting the expected rare semileptonic background. The 
interpolation procedure assumes a flat background shape and has a systematic uncertainty of 
4%, which is evaluated by varying the flight-length significance selections and by using a linear 
background shape with a variable slope. 
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Figure 5: B+ — > J/t/'K+ invariant-mass distributions in the barrel (left) and endcap (right) chan- 
nels. The solid (dashed) lines show the fits to the data (background). 

Figure |6] shows the measured dimuon invariant-mass distributions. In the sidebands the ob- 
served number of events is equal to six (seven) for the barrel (endcap) channel. Six events 
are observed in the Bg — ?► }i^}i^ signal windows (two in the barrel and four in the endcap), 
while two events are observed in the B*^ — ?► }i^}i^ barrel channel and none in the endcap chan- 
nel. As indicated by the numbers shown in Table [2| this observation is consistent with the SM 
expectation for signal plus background. 



Upper limits on the Bg — ?► fi^}i^ and B^ — ?► fi^fi^ branching fractions are determined using the 
CLg method [21 , 22J. Table|2]lists all the values needed for the extraction of the results for both 
the barrel and endcap channels. The combined upper limits for the barrel and endcap channels 
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^-9\ 



are B(B^ -^ H^H') < 7.7 x 10"^ (6.4 x 10"^) and B(B" -^ ^+M") < 1-8 x 10"^ (1.4 x 10" 
at 95% (90%) CL. The median expected upper limits at 95% CL are 8.4 x 10"'^ (1.6 x 10"'^) for 
Bg — > }i^}i^{B^ — > }i^}i^), where the number of expected signal events is based on the SM 
value. Including cross-feed between the B^and Bg decays, the background-only p value is 0.11 
(0.24) for Bg — > }i^}i^{?P —> ji+ji^), corresponding to 1.2 (0.7) standard deviations. The p value 
for the background plus SM signal hypotheses is 0.71 (0.86) for Bg 



^+^-{B^ 



¥^V )■ 



Table 2: The event selection efficiency for signal events £tot/ the SM-predicted number of signal 
events N^^^ai' ^^ expected number of peaking background events Np^^j, and combinatorial 
background events N^^^^^z and the number of observed events Nobs in the barrel and endcap 
channels for Bg — ?> fi^ji^ and B*^ — > }i^H^ ■ The quoted errors include both, the statistical and 
the systematic uncertainties. 



Variable 


B'J ^ ^+jM- Barrel 


B^ -^ ]A+li- Barrel 


B^ ^ li+]A- Endcap 


B^ -^ /^+/^- Endcap 


£tot 


0.0029 ± 0.0002 


0.0029 ± 0.0002 


0.0016 ± 0.0002 


0.0016 ± 0.0002 


n;Li 

'comb 


0.24 ± 0.02 
0.33 ± 0.07 
0.40 ± 0.34 


2.70 ± 0.41 
0.18 ±0.06 
0.59 ± 0.50 


0.10 ±0.01 
0.15 ±0.03 
0.76 ± 0.35 


1.23 ±0.18 
0.08 ± 0.02 
1.14 ±0.53 


-"^total 


0.97 ±0.35 


3.47 ± 0.65 


1.01 ±0.35 


2.45 ± 0.56 


Nobs 


2 


2 





4 



H^¥' 



and B*^ — > ji^fi has been performed 



7 Summary 

An analysis searching for the rare decays Bg 
in pp collisions at ^/s = 7 TeV. A data sample corresponding to an integrated luminosity of 
5fb^ has been used. This result supersedes our previous measurement [9J. Stricter selection 
requirements were applied, resulting in a better sensitivity and a higher expected signal-to- 
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background ratio. The observed number of events is consistent with background plus SM 
signals. The resulting upper limits on the branching fractions are S(Bg — )■ fi^fi^) < 7.7 x 10^^ 
and 5(5*^ -^ /^^f') < 1.8 x 10^^ at 95% CL. These upper limits can be used to improve bounds 
on the parameter space for a number of potential extensions to the standard model. 
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